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The low energy phonons of two different graphite intercalation compounds (GICs) have been 
measured as a function of temperature using inelastic x-ray scattering (IXS). In the case of the 
non-superconductor BaCe, the phonons observed are significantly higher (up to 20%) in energy 
than those predicted by theory, in contrast to the reasonable agreement found in superconducting 
CaCe. Additional IXS intensity is observed below 15 meV in both BaCe and CaCe- It has been 
previously suggested that this additional inelastic intensity may arise from defect or vacancy modes 
not predicted by theory (d'Astuto et al, Phys. Rev. B 81 104519 (2010)). Here it is shown that 
this additional intensity can arise directly from the polycrystallinc nature of the available samples. 
Our results show that future theoretical work is required to understand the relationship between 
the crystal structure, the phonons and the superconductivity in GICs. 
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I. INTRODUCTION 

Since the discovery of superconductivity in YbCe (Tg = 
6.5 K) and CaCg {T^ = 11.5 K) at temperatures over 
an order of magnitude higher than previously found in 
graphite intercalation compounds (GICs)!'^, the proper- 
ties of this family of GICs have been extensively stud- 
ied by a variety of different experimental techniques'^"^. 
Although initially an exotic superconducting mechanism 
was proposed involving acoustic plasmons^, subsequent 
density functional theory (DFT) studies described the 
superconductivity via a more orthodox electron-phonon 
(e-ph) coupling mechanism with s-wave symmctry^""'^^. 
These DFT descriptions predict that the e-ph coupling 
is approximately equal for phonons associated with vi- 
bration of the carbon atoms and for phonons associated 
with movement of the intercalant calcium. 

At present there are experimental studies in the lit- 
erature which give conflicting viewpoints concerning the 
nature of the e-ph coupling in GICs. A large Ca iso- 
tope effect (a(Ca) ~ 0.5) has been measured in CaCe^, 
which, if viewed within the BCS description of super- 
conductivity, suggests that only the phonons due to the 
vibration of calcium are involved in the electron pairing. 
In contrast, angle-resolved photoemission spectroscopy 
(ARPES) measurements on CaCe* have found that the 
e-ph coupling to graphite-like high-enc!rgy phonons is so 
strong that it can explain the superconducting transition 
temperature alone, without any additional coupling to 
calcium phonons. These discrepancies point to the need 
for a detailed study of the phonons in GICs, to both 



test the DFT description and to look for direct evidence 
for e-ph coupling involving specific phonons. Moreover, 
phonon studies in graphitic systems in general are impor- 
tant as the electron-phonon interactions in these systems 

are under much scrutiny "'^^""'^®. 

Superconductivity in GICs is directly linked to the 
graphite layer separation, with the superconducting tran- 
sition temperature Tc increasing as the graphite layer 
separation d is reduced. This trend is supported both 
by the observed values of Tc in a variety of GICs''' and by 
measurements of Tc as a function of pressure in CaCe^^ 
and YbCe^*. Indeed in BaCe the graphite layer sepa- 
ration is so large that superconductivity appears to be 
suppressed entirely, with no superconducting transition 
observed down to 0.080 K^^. It is therefore instructive to 
study the STiperconductivity in GICs by studying GICs 
with difl'erent intercalants, since by changing the inter- 
calant one changes d and therefore tunes Tc. 

In this paper we present the low energy phonon disper- 
sions in non-superconducting BaCe and superconduct- 
ing CaCe as measured using inelastic x-ray scattering 
(IXS) . These data represent the first momentum-resolved 
phonon measurements on BaCe. We find a substan- 
tial discrepancy between experiment and theory in the 
phonon energies of BaCe, in contrast to the good agree- 
ment in the case of CaCe. Like many other layered ma- 
terials, GICs are difficult to synthesise as high-quality 
single crystals. We demonstrate here that the details of 
the preferred orientation (texture) of the crystallites in 
these GIC samples can lead to the observation of phonons 
which, because of their polarization, are theoretically for- 
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bidden to be observed. Our work highlights the impor- 
tance of accounting for these effects in phonon studies of 
textured polycrystalhne samples and calls for the need 
for single crystal GIC samples to address the role of the 
electron-phonon coupling in these materials. 

The low energy phonons in CaCg have been measured 
previously^" and found to be in good overall agreement 
with the published DFT calculations. Subsequent to this 
work two momentum-resolved phonon studies have been 
made on CaCe, an inelastic neutron scattering (INS) 
study, which concentrates primarily on the high energy 
graphite-like phonon modes^^, and a study of the low 
energy phonons in CaCe performed using both IXS and 
INS^^. While in Ref.^^ good agreement was found with 
the calculated DFT phonon dispersions by taking ac- 
count of the polycrystalhne nature of the sample, in 
Ref.^^ it was suggested that an additional phonon mode 
exists in CaCg of uncertain origin. Here we propose an 
alternative explanation: that the additional IXS inten- 
sity arises from the weak crystallographic texture in the 
polycrystalhne GIC samples. 
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II. EXPERIMENTAL METHODS 

The GIC samples were made using ZYA grade highly 
orientated pyrolytic graphite (HOPG) platelets pur- 
chased from GE Advanced Ceramics. The BaGg sam- 
ple was made using the vapour transport method^'^. An 
HOPG platelet was outgassed at 500° C and then sealed 
into a quartz tube along with the barium metal (purity 
> 99.99%) under high vacuum (< 10~^ mbar). The tube 
was then heated to 490° C and maintained at this tem- 
perature for 4 weeks. The CaCg samples used were pre- 
pared by immersing a HOPG platelet in a Li-Ca alloy for 
10 days, as described elsewhere^''. GaCg samples made 
from the same batch were found to have a sharp super- 
conducting transition at 11.5 K from magnetic suscepti- 
bility measurements. The very high purity of the sam- 
ples (> 99% pure in both cases) can be seen in the {001) 
diffraction shown in Figure 1 (e) and 1 (f), where there 
are no visible Bragg peaks from any impurities. 

The starting graphite (HOPG) is composed of small 
crystallites (~ 1 /im) which have a strong preferred ori- 
entation (strong texture) perpendicular to the graphene 
planes (out-of-plane) , giving a (OOZ) mosaic with full- 
width half-maximum (FWHM) as low as 0.2°. However 
within the graphene planes (in-plane) , the crystallites are 
orientated randomly^'^. After intercalation the crystal- 
lites are still oriented randomly in-plane, but out-of-plane 
the orientation of the crystallites is more random, with 
the GICs studied here having (001) mosaics of 5°. This 
means that the texture is weaker in the GIC samples than 
in HOPG, since a weaker texture means that the samples 
are more like a perfect powders, which have zero texture. 

CaCe has the structure R3m, which can be described 
using a rhombohedral or a hexagonal basis^. The cal- 
cium atoms are arranged in three different ways in dif- 



FIG. 1; The structural properties of BaCe and CaCe: (a)-(b) 
Crystal structure, (c)-(d) First Brillouin zone with symmetry 
points, (e)-(f) Diffraction on BaCe and CaCe measured in 
{001) direction at 5 K (dashed blue line) and at 300 K (solid 
black line) at 3-ID. The symmetry points shown are at the 
edge of the Brillouin zone apart from the F point, which is at 
the zone centre. There is practically no signal from impurities 
in the {001) diffraction from either sample. 



ferent intercalant layers (called AaA/3A7 stacking, where 
the Roman letters define graphite layers and the Greek 
letters intercalant layers). The unit cell of CaCg in the 
hexagonal basis is shown in Figure 1(b), and the shape of 
the first Brillouin zone is shown in Figure 1(d). Here we 
define the reciprocal lattice directions using the hexago- 
nal basis, meaning that the out-of-plane direction is the 
{001) direction. This convention aids comparison with 
the BaCg data, as BaCg has the space group Pd^/mmc, 
which is normally described within the hexagonal basis. 
The stacking in BaCe is AaA/3, as shown in Figure 1(a). 
The first Brillouin zone of BaCe is presented in Figure 
1(c). The lattice parameters of CaCg are a = 4.333(2) A 
and c = 13.572(2) A^^, giving a graphite layer separation 
d of 4.524(1) A, and in BaCe the lattice parameters are 
a = 4.302(6) A and c = 10.50(4) A with d = 5.25(2) A^^. 

The IXS measurements were performed at 3-ID at the 
Advanced Photon Source, Argonne National Laboratory, 
with an incident x-ray energy of 21.657 keV^^, providing 
an energy resolution between 2.2 and 2.4 meV, depending 
on the specific analyzer. The spectrometer has four ana- 
lyzers, allowing energy scans to be done at different mo- 
mentum transfers simultaneously. The momentum reso- 
lution in all cases was 0.072 A~^ in the scattering plane 
and 0.183 A~^ perpendicular to it. The phonon peaks 
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were fitted with pseudo-Voigt functions, which were ap- 
propriately scaled by the Bose factor. The phonon dis- 
persions were produced by plotting the peak positions as 
a function of the phonon wavevector q, defined within the 
equation Q = G + q, where G is the nearest reciprocal 
lattice vector and Q the momentum transfer. 

In order to model the effect of the crystallographic tex- 
ture in the GIC samples on the phonon spectra, simula- 
tions of the IXS data were produced by summing hun- 
dreds of simulated IXS intensities, each of which was per- 
formed at a specific momentum transfer. The crystallo- 
graphic texture was described by performing simulations 
over a volume in reciprocal space expressed in spherical 
polar coordinates {\Q\,9,(j)), where 9 has its rotation axis 
out-of-plane. For the 2D powder simulations, the mo- 
mentum transfers were selected using Lorentzian sam- 
pling of (p with a FWHM of 5° and allowing 9 to take 
any value. For the 3D powder simulations both and 9 
were allowed to take any value. The summed IXS spectra 
were then convolved with the momentum and energy res- 
olution of the IXS spectrometer. This method was also 
used in our recent INS study of CaCe^^. 

The two-dimensionality of these GICs, together with 
the significant difference in mass between the intercalant 
and carbon, means that the phonon modes can be sep- 
arated, to a good approximation, into four groups: Ixy, 
Iz I and Cz , where l^y describes phonon modes purely 
due to vibrations of intercalant atoms in-plane, I^ the in- 
tercalant phonons out-of-plane, and Cxy{z) the equivalent 
carbon in-plane (out-of-plane) phonons. 
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FIG. 2: BaCe (001) phonon dispersion measured at 5 K 
(squares) and 300 K (circles). Empty symbols denote Ba^ 
phonon intensity, full symbols label Ba^y phonon intensity. 
The Ba^jy phonons should not be visible in this scattering ge- 
ometry (see text). The theoretical dispersion of the Baz and 
Ba,xy phonons are plotted with solid lines and dashed lines 
respectively^^. A guide to the eye is plotted showing the dis- 
persion of the Baz phonons at 300 K (thick solid line). The 
symmetry point A is located at g = 0.299 A~^. The crosses 
mark the phonon energies calculated at F using the experi- 
mental lattice parameters and space group of BaCe . 



III. INELASTIC X-RAY SCATTERING 
MEASUREMENTS 

A. Phonons in BaCs 

The (001) phonon dispersions of BaCg measured at 
both 5 K and 300 K are plotted in Figure 2, together 
with calculated dispersions^^. The measured Baz phonon 
branches are significantly higher in energy than the theo- 
retical values: near the edge of the Brillouin zone, at sym- 
metry point A, the measured phonon energy for the optic 
Ba^ mode is almost 20% larger than that predicted. In 
addition, in all of the measured IXS energy scans on BaCg 
there is additional intensity observed below 10 meV, plot- 
ted with full symbols, which is not predicted by theory. 

So what is the origin of this discrepancy between ex- 
periment and theory in BaCg? We consider a number of 
possibilities. The DFT calculations were performed us- 
ing a different structure to the structure experimentally 
determined, which may have affected their results. In 
addition, the charge transfer from the Ba atoms to the 
graphene planes may be inaccurately predicted by theory. 
Finally, the calculations do not account for the polycrys- 
tallinity of the real samples, so part of the disagreement 
may be due to their crystallographic texture. 

The published calculations for BaCg use lattice param- 



eters equivalent to a = 4.350 A and c = 10.40 A^^ which 
are significantly different to the values found via x-ray 
diffraction (a = 4.302(6) A and c = 10^50(4) A)^^. In 
addition, the space group of CaCg (i?3m) is used for 
BaCg in the calculation, rather than the experimentally 
found PQ^/mmc. To determine whether the observed 
discrepancy arises because of these structural differences, 
we performed an additional phonon calculation at the F 
point using the experimental structure and lattice pa- 
rameters, as shown in Figure 2 (crosses). This calcu- 
lation gives phonon energies in approximate agreement 
with the previous calculation, showing that the calcula- 
tions are largely insensitive to small changes in both the 
lattice parameters and the space group. Therefore the 
incorrect structure used in the initial calculations can be 
eliminated as a cause of the discrepancies. 

Another possibility is that the charge transfer from 
the Ba atoms to the graphite has been underestimated 
theoretically. If there is less charge than predicted in the 
graphitic tt* bands, then the bonds should be stronger 
than predicted, since filling the anti-bonding tt* band 
destablizes the bonds. Stronger bonds lead directly to 
higher phonon energies. The effect of charge transfer on 
phonons in GICs has been predicted theoretically^^ and 
observed using Raman scattering^^. This would explain 
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why the phonons are higher in energy than predicted, 
but cannot explain why additional intensity is observed 
below 10 meV. 

The additional phonon intensity can be explained as a 
result of the crystallographic texture of the GIG samples. 
Here we argue that this signal arises from the excitation 
of Ixy phonons over a large volume in reciprocal space 
due to the weak crystallographic texture. For an ideal 
single crystal, and perfect instrumental resolution, the 
Ixy phonons are disallowed (their intensity is zero) when 
Q is completely out-of-plane due to the Q ■ e{q) term in 
the IXS phonon cross-section, where e{q) is the eigen- 
vector of the phonon with wavevector q^^. However in 
the case of the polycrystalline GIG samples studied, the 
weak crystallographic texture provides an explanation for 
the observation of the l^y phonons. Even if the nominal 
momentum Q is entirely out-of-plane, the large recipro- 
cal space volume integrated over in each measurement 
due to the weak preferred orientation of the crystallites 
will include many values of Q which have a significant 
component in-plane. This effect is discussed further in 
Section IV. 



B. Phonons in CaCe 

Figure 3 presents the (001) phonon dispersion in GaGe 
at 5 K, 50 K and 300 K. A subset of the GaGg data 
has already been published^", but the scope of the data 
presented here is much more extensive. The additional 
phonon intensity below 15 meV again results from the 
large mosaic of the sample and is discussed at length 
in Section IV. The phonon dispersions calculated using 
DFT are plotted on the same figure^°. The two Gaz 
modes are well described by the DFT calculations over 
the whole Q range sampled, especially in the case of the 
acoustic Gaz mode. The higher energy mode which dis- 
perses between 30 and 40 meV (the optic Gaz mode) is 
about 2 meV higher in energy than predicted, but the 
character of the dispersion is reasonably well described. 

In both GaGg and BaGg the energies of the Iz modes 
are slightly hardened (< 1 meV) upon cooling from 300 
K to 5 K, but there is no observable difference in GaGg 
between the data measured above and below Tc- The 
small temperature dependence most likely results from 
the reduction in the c lattice parameter upon cooling, 
visible in the diffraction presented in Figure 1(e) and 1(f). 



IV. MODELLING THE CRYSTALLOGRAPHIC 
TEXTURE IN BaCe AND CaCg 

In Figures 4 and 5 a selection of phonon spectra mea- 
sured at 300 K in the (001) direction in BaGg and GaGg 
are presented. In both cases two phonons are observed 
in the raw IXS data (the acoustic and optic Iz branches) , 
as well as additional IXS intensity at low energies. In 
each of the panels two IXS simulations are plotted. The 
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FIG. 3: CaCe (001) phonon dispersion measured at 5 K 
(squares), 50 K (triangles) and 300 K (circles). Empty sym- 
bols denote Ca^ phonon intensity, full symbols label Ca^y 
phonon intensity. The C&^cy phonons should not be visible in 
this scattering geometry (see text). The theoretical disper- 
sion of the Ca.z and Ca^y phonons are plotted with solid lines 
and dashed lines respectively^". A guide to the eye is plotted 
showing the dispersion of the experimental Gaz phonons at 
300 K (thick solid line). The symmetry point L is located at 
q = 0.694 A-\ 



first simulation models the sample as a 2D powder: that 
is, the crystallites have no preferred orientation in-plane, 
but have a Lorentzian mosaic of FWHM — 5° out-of- 
plane, consistent with our x-ray diffraction. The second 
simulation models the sample with no preferred orienta- 
tion (no texture): i.e. as a 3D powder. Both simulations 
are added to the experimental elastic intensity. 

In the case of BaGg, plotted in Figure 4, the overall 
agreement is better with the 2D powder model than with 
the 3D powder model, especially in the ratio between the 
acoustic and optic Baz phonons. However in the 2D pow- 
der model the IXS intensity due to the Baa;^ phonons is 
much smaller than measured. The additional features 
predicted by the 3D powder model are similar to those 
observed, but there the intensities of the Baz phonons 
are underestimated. These observations suggest that per- 
haps the crystallographic texture of the BaGg samples is 
weaker than expected from our x-ray diffraction, lying 
somewhere between the 2D and 3D powder models. 

In the GaGg data in Figure 5, the low energy IXS fea- 
tures look qualitatively similar to the CsL^y features pre- 
dicted by the 3D powder model. However the 3D pow- 
der description of GaGg is once again not a satisfactory 
description, as the intensity of the Gaz phonons is un- 
derestimated in most cases and the IXS intensity at low 
energies is significantly overestimated'^'^. The measured 
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FIG. 4: (a) (00/) IXS data at g = 2.91 A"^ on BaCe at 
300 K, plotted with the data fit and the decomposed fitted 
phonon peaks, (b) Fit to the data at Q = 2.91 A~^, plotted 
with the fitted elastic line and against IXS simulations based 
on ab-initio calculations^^ which use the 2D powder model 
and the 3D powder model described in the text, (c) & (d) 
Identical to (a) & (b) but for Q — 3.12 A~^. In each panel the 
experimental (simulated) IXS intensity due to Ba^j, phonons 
is marked with black (blue) arrows. 



CaCe data appear to lie somewhere between the 2D and 
3D powder simulations, similar to the BaCe data, though 
the CaCe appears to be more like the 3D powder than 
BaCg. This suggests that the distribution of intercalant 
atoms in CaCg is rather random, which is consistent with 
the behavior of CaCg under pressure as studied with 
x-ray diffraction, where the Ca atoms are found to be 
very mobile^. More theoretical and experimental work 
is required in order to understand the complicated crys- 
tallographic texture of GICs, with more extensive x-ray 
diffraction being a natural starting point. 

A recent IXS and INS study by d'Astuto et al.^^ on 
CaCe suggested that the additional inelastic intensity 
was due to an interaction with the acoustic Ca^ mode, 
causing an avoided crossing, or anti-crossing, which is 
seen as a splitting in the acoustic Ca^ mode. Their work 
was supplemented by INS data, which allowed them to 
more easily access momenta nearer to the T point. The 
study concluded that the additional inelastic intensity 
could be due to a defect or vacancy mode. Although we 
cannot exclude this hypothesis, our simulations suggest 
that if the orientation of the crystallites in these sam- 
ples is more random than previously thought, the weak 
crystallographic texture can account for the anomalous 
features without recourse to such a mode. 



V. CONCLUSIONS 
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FIG. 5: (a) (00/) IXS data at Q = 2.42 A~i on CaCe at 
300 K, plotted with the data fit and the decomposed fitted 
phonon peaks, (b) Fit to the data at Q = 2.42 A~^, plotted 
with the fitted elastic line and against IXS simulations based 
on ab-initio calculations^" which use the 2D powder model 
and the 3D powder model described in the text, (c) & (d) 
Identical to (a) & (b) but for Q — 2.08 A~^. In each panel the 
experimental (simulated) IXS intensity due to Caa;j, phonons 
is marked with black (blue) arrows. 



To summarise, the dispersions of the low-energy 
phonons in BaCg and CaCg have been measured as a 
function of temperature using inelastic x-ray scattering. 
In BaCe the experimental and DFT-calculated phonon 
dispersions^^ disagree, with measured phonon energies 
up to 20% higher than predicted. We suggest that this 
large discrepancy may result from an underestimation 
of the charge transfer from the Ba atoms to the graphite 
sheets in the theory. Our work motivates further study on 
BaCe in order to examine the underlying reasons for this 
disagreement. In contrast, reasonable agreement with 
theory is found in CaCg for the Ca^ phonons. This con- 
sistency between theory and experiment provides indirect 
supporting evidence for the DFT description^'' of the su- 
perconductivity in CaCe- 

No signatures of electron-phonon coupling are observed 
in the phonon dispersions or the phonon widths in either 
non-superconducting BaCe or superconducting CaCe, 
despite the drastically different superconducting transi- 
tions of these related compounds. In both BaCg and 
CaCg there is a small (< 1 meV) hardening of the Iz 
phonons as the temperature is decreased, but this is likely 
due to a reduction in the c lattice parameter and is un- 
affected by the presence of superconductivity in CaCg 
below 11.5 K. The largest source of phonon broadening 
experimentally is very likely the weak crystallographic 
texture inherent in the GIG samples. 

Finally, the IXS simulations presented here show that 
weak crystallographic texture in polycrystalline GIG 
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samples may lead to additional inelastic intensity from 
Ixy phonons. Such additional intensity has been observed 
recently in CaCg^"'^^ and YhCe^^, but also in older INS 
studies on RbC24'^^ and KC24^^. This work provides a 
timely reminder that the crystallographic texture inher- 
ent in many graphitic systems may give rise to unex- 
pected experimental effects. 
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